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Abstract 

The LIT approach is reviewed both for inclusive and exclusive reactions. It is shown that the 
method reduces a continuum state problem to a bound-state-like problem, which then can be 
solved with typical bound-state techniques. The LIT approach opens up the possibility to perform 
ah initio calculations of reactions also for those particle systems which presently are out of reach 
in conventional approaches with explicit calculations of many-body continuum wave functions. 
Various LIT applications are discussed ranging from particle systems with two nucleons up to 
particle systems with seven nucleons. 
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I. INTRODUCTION 



Ab initio calculations are a central element of few-body physics. The only input in such 
calculations is a well defined Hamiltonian. In nonrelativistic nuclear physics one then solves 
the Schrodinger - or equivalent - equations without introducing any approximation. In 
order to calculate cross sections for reactions, where final or initial states are scattering 
states in the continuum, one is faced with the problem of an ab initio calculation of such 
continuum states. It is well known that already the calculation of a three-body continuum 
wave function is quite difficult and that today a complete four-body calculation, with all 
possible break-up channels open, is out of reach. However, the problem of a many-body 
continuum state calculation can be circumvented if one uses the Lorentz integral transform 
(LIT) method In fact the LIT approach allows the ab initio calculation of reaction 
cross sections, where a many-body continuum is involved, without requiring the knowledge 
of the generally complicated many-body continuum wave function. The scattering problem 
is reduced to a calculation of a localized function with an asymptotic boundary condition 
similar to a bound-state wave function. Such an approach was already proposed by Efros 
in 1985, but with the Stieltjes instead of the Lorentz integral transform 2|. However, it has 
been found that the application of the Stieltjes transform is problematic since it leads to 
serious inversion problems 3|. 

The LIT method has been applied to various electroweak cross sections in the nuclear 
mass range from A=3 to A=7. Among the applications are the first realistic ab initio cal- 



culations of the nuclear three- and four-body total photoabsorption cross sections 
as well as of the inelastic neutral current neutrino scattering off ^He jol. In addition first 
ab initio calculations have been performed for the total photoabsorption cross sections of 
^'^He and ^'^Li with semirealistic forces [?, y, Q]. Other applications were carried out for 



the inelastic inclusive electron scattering cross section (see e.g. 1^, lUl, ll^ ) . Besides inclu- 



sive electroweak reactions also LIT calculations of exclusive reactions have been performed, 
namely for ^Re{'y,n)^Re and ^Re{'y,pfR Q, ^Re{e,e'pfR Q, and ^Re{e,e'd)d Q- Fur- 
ther applications and a detailed description of the LIT method are presented in a recent 
review article 16 1. 



The paper is organized as follows. In section II the general form of inclusive and exclusive 
cross sections of a particle system induced by an external probe is outlined. Sections HI and 
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IV describe the LIT approach for inclusive and exclusive response functions, respectively. 
In section V various LIT applications are discussed. 



II. STRUCTURE OF ELECTROWEAK CROSS SECTIONS 

Perturbation-induced reactions can be divided in inclusive and exclusive processes. In 
the former case the final state of the particle system is not observed. In the latter case the 
final state is at least partially observed, like e.g. in an {e,e'p) reaction, where besides the 
scattered electron also an outgoing proton with energy Ep and scattering angle flp = {9p, $p) 
is detected. 

Inclusive cross sections of perturbation-induced reactions have the following general form 



(see e.g. [ll 

M 



ttext 



'^fi{uj,q,eext)Ri{uJ,q), (1) 



dudQext . , 

1=1 

where u and q and are energy and momentum transfer of the external probe to the particle 
system, Qext = {^ext, (pext) denotes the scattering angle of the external probe, Ocxt is a 
constant characteristic for the external probe, and fi are kinematic functions. The functions 
Ri describe the various responses of the particle system to the external probe and thus 
contain information about the dynamics of the particle system. They are defined as follows 



R.{u,q) = ^\{f\em'Sico-iEf-Eo)). 



(2) 



Here Eq and |0) are ground state energy and wave function of the particle system under 
consideration, Ef and |/) denote final state energy and wave function of the final particle 
system, and 6j is the operator inducing the response function R^. 

As an example for an exclusive cross section we consider the (e, e'p) case, here one has 

d^a ^ 

dwdn dn " ^ ^' ^e)^i(0p)^i(^> Q, Ep) , (3) 

^ P i=l 

where the (^p dependence of the cross sections is described by the known functions gi{(t>p). 
Exclusive response functions do not have such a simple form as the inclusive functions Ri. 
For their definition we refer to here we only mention that transition matrix elements 
from the ground state \0a) to a specific final state 

T-f^, = {fm\Oa), (4) 
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are their essential ingredients, where a and (5 stand for additional quantum numbers of 
initial and final state wave functions. 

The following relation between the inclusive Ri of ([T]) and the exclusive rj of ([3]) holds: 

Q) = j dnpdEpgi{(j)p)ri{u;, q, 9p, Ep) . (5) 

Note that the number of exclusive response functions is generally greater than the number 
of inclusive ones (A^ > M), since the integration over the azimuthal angle of the outgoing 
particle can yield zero, i.e. the integration over 0p in the {e,e'p) example above. 

With additional polarization degrees of freedom for beam and/or target and/or outgoing 
particles many more additional inclusive and exclusive response functions can be defined 
(see e.g. the deuteron case in [l8|). 



III. CALCULATION OF INCLUSIVE RESPONSES WITH THE LIT METHOD 

As already mentioned, with the LIT method one avoids the explicit calculation of scatter- 
ing wave functions. Instead, for the calculation of Ri of ([2]) one proceeds in the following way. 
One first calculates the ground state wave function |0) of the particle system in question. 
Then one solves the equation 

{H-Eo-aR-tai)\^,) = e,\0), (6) 

where H is the Hamiltonian of the particle system and a^/j are parameters, whose meaning 
is explained below. Since the eigenvalues of H have to be real the homogeneous version of 
([6]) has only the trivial solution = and thus ([6]) has a unique solution. In addition, due 
to the asymptotically vanishing ground state wave function |0) also the right-hand-side of ([6]) 
vanishes asymptotically. Therefore, and because of the complex energy Eq — a^ — iaj, has 
a similar asymptotic boundary condition as a bound state. It means that is a so-called 
localized function, i.e. square- integrable with a norm This has very important 

consequences: even if the aim is a calculation of a reaction cross section in the continuum, 
one is not confronted with a scattering state problem any more, in fact one needs to apply 
only bound-state techniques for the solution of ([6]). 

The key point of the LIT method consists in the fact that the Lorentz integral transform 
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Li of the response function 

U{o-R, ai, q) = J Rii^^ q) (^R^ CT/) (7) 

is related to the norm which can be obtained from the solution of In fact one 

has 

Li((Ji?,cr/,g) = {^i{aR,ai,q)\^i{(TR,(Ti,q)) (8) 

(g dependence of Lj and \E'j will be dropped in the following). In ([7]) £ is a Lorentzian 
centered at ctr with a width F = 2a i: 

C{uj, aR, crj) = ^— — 2 . (9) 

{iu - aRY + aj 

Now also the meaning of the parameters o"/?// becomes evident: aj represents a kind of 
energy resolution, while with aR a, given energy range can be scanned. 

With the above equations the principle idea of the LIT method can be explained: one 
solves the LIT equation ([6]) for many values of aR and a fixed o"/, calculates hi{aR, aj =const), 
and then one inverts the transform in order to determine Ri{uJ, q). 

Before coming to the inversion we first want to derive the relation ([8]). Starting from the 
definition of Lj in ([7]) one has 

J [u - aRY + aj J [uj- aR + ai) [u - aR- aj) 
Using ([2]) and carrying out the integration in du one gets 

, 2:;(o|ell/)(/|9.|o) 

duj- — -5{uj - {Ef - Eo)) 

{uJ -aR + aj) {u - aR - aj) 



t 



m\{Ef -Eo-aR + taj)-'\f){f\{Ef -E^-aR- iaiY^Q,\Q). (11) 
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Then one replaces Ef hy the Hamilton operator H and uses the closure property of the 
eigenstates of the Hamiltonian I/) (/I = 1): 

L, = {Q\Q\[H~E^~aR + laIY\H-E^-aR-^aY-'Q^\^) 

= {^^\^^) (12) 

with 

l^j) = {H-Eo-aR- ^aIY'&^\0) • (13) 
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One sees that relation ([8]) is indeed obtained and that fulfills (|6]). 

The standard LIT inversion method consists in the following ansatz for the response 
function 

Mm ax 

Ri{uj\q) = ^ CmXmi^^' , aj) , (14) 



m=l 



here the argument u of Ri is replaced hj u' = uj — ujth, where Uth is the break-up threshold 
of the reaction into the continuum. In case of LIT contributions attributed to bound states 
due to, e.g., elastic transitions, one can easily subtract such contributions in order to obtain 



an "inelastic" LIT (see 



16|). The Xm are given functions with nonlinear parameters aj 



Normally the following basis set is taken 



xUuj',a,)=u;'^^exp{-^). (15) 



Substituting such an expansion into the right hand side of ([7]) one obtains 



M, 



max 



Ui(^R,cri) = ^ CmXm(cri?,o-/,«i) , (16) 

m=l 

where 

Xm{<yR,(^i,aj) = I du' -^"^^^'"^^ . (17) 
Jq [uj' - anY + cr| 

For given values of aj and Mmax the linear parameters Cm are determined from a best fit 
of hi{aR,ai) of (fT6!) to the calculated Lj((TR,cr/) of ([8]) for a fixed o"/ and a number of 
points much larger than Mmax- In addition one should vary the various nonlinear parameter 
Uj over a sufficiently large range. The parameter ai, however, can in general be determined 
from the known threshold behavior of the response function. One has to increase M^ax up 
to the point that a stable inversion result is found for some range of Mmax values, which 
then can be taken as final inversion result. Note, however, that a too large value of Mmax 
might lead to an oscillatory behavior of R^. The origin for such an unphysical behavior lies 
in the precision of the calculated Lj. If the precision is further increased, higher and higher 



^max values can in principle be used in the inversion (see also 19|). 

One can repeat the whole procedure with a second o"/ value. Of course, the inversion 
should lead to the same Ri result as with the previous ct/. The basis set Xm can also be 
modified in order to take into account narrow structures like resonances (see section V.A). 
More information concerning the inversion and alternative inversion methods are found in 

3- 



IV. CALCULATION OF EXCLUSIVE RESPONSES WITH THE LIT METHOD 



For the exclusive response function rj one has to evaluate T-matrix elements of the type 
given in (jl]). One starts the LIT calculation using the general form of the final state wave 



function for the considered break-up channel 20|, 



\^j{Ej)) = \^f{Ef)) + {Ef-H- ii^)-^V\^f{Ef)) , (18) 

where \^{Ef)) is a so-called channel function (with proper antisymmetrization) given in 
general by the fragment bound states times their relative free motion and V is the sum of 
potentials acting between particles belonging to different fragments. Thus the transition 
matrix element Tq/^j (additional quantum numbers are dropped) takes the following form 

To;, = (vl/(E/)|e.|0) 

= mEf)\e,\0) + mEf)\V{Ej-H + tr^)-'\0). (19) 

The first term of the right-hand-side is the so-called Born term (TJy-°™), which can be evalu- 
ated without greater problems. The second term [T^^j) depends on the final state interaction 
and its evaluation is much more difficult. However, using the LIT approach, one can proceed 
as follows. One rewrites T^^j in a spectral representation, 

/-oo 

= J2(Ej - EMM, En) + {Ej-E' + ir^r^FoUEf, E')dE' (20) 

with 

FofAEf, E') = ^ {^{Ef)\V\m,){m,\Q,\mE} - E') . (21) 

The function F^j^i has a similar form as an inclusive response function therefore one can 
apply an analogous LIT method as in the inclusive case, however left- and right-hand sides 
are not identical, hence two LIT equations are obtained: 

[H-an- 2a/)|^,) = e,|0) , {H-aR- iai)\^>y) = V\^{Ef)) . (22) 

The first one is essentially the same as ([6]). The second equation has a different right-hand 
side, but with the important feature to vanish asymptotically for a nuclear potential V . 
Therefore the equation can again be solved with a bound-state technique. In case of an 
additional Coulomb interaction, one may use Coulomb wave functions instead of the free 
motion if only two of the fragments carry charge. 



Having calculated and one evaluates the overlap (\Efy|^i); which is identical to the 
LIT of the function -fo/,*; and hence Fo/_j is obtained by the inversion of the LIT. The FSI 
part of the T-matrix element is then given by 

POD 

T^SiEf) = -^vrFo/,(£;/, Ej)+V {Ej - E')-'FoUEj. E')dE' (23) 

and the sum of T^j'] and the simpler Born term T(^°™ leads to the total result for the 
transition matrix element. 

As shown in [16i] the LIT formalism for exclusive reactions can be reformulated such that 
only a solution for is requested, while is not needed. Both possibilities have been 
used in 2l(|, where the d{e, e'p)n reaction has been calculated as a test case for the exclusive 
LIT formalism. 



V. APPLICATION OF THE LIT METHOD 

As we have shown in the previous sections the main point of the LIT approach consists 
in the fact that a scattering state problem is reduced to a bound-state-like problem. In 
other words the calculation of continuum wave functions is not required, instead one has 
to solve equations which can be solved with bound-state techniques. For A>2 the calcula- 
tion of continuum wave functions is difficult or today even impossible, thus, with the LIT 
method, one can extend the range of calculations to considerably larger A. In fact one may 
conclude the following: if one is able to carry out a bound-state calculation for a given 
particle system then the LIT approach opens up the possibility to perform calculations for 
continuum reactions with this particle system. In principle one is not restricted to use a 
specific bound-state technique, but in most LIT calculations an expansion of ground state 
|0) and LIT function |\E') in hyperspherical harmonics (HH) is employed. Information con- 



cerning such expansions is given in 



16j, here we only want to mention that the realistic 



LIT applications for A=3 have been performed with the CHH technique 



221 ]. whereas the 



realistic (semirealistic) applications for A=4 (A>4) have been carried out with the EIHH 



approach [23 1. 



For the solution of the LIT equation (I6l) the Lanczos method is used in most cases [24 1. 
In this context it should be pointed out that the LIT method is different from an approach 
where a so-called Lanczos response -RLanc is introduced, which is essentially a LIT with small 
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0"/, which, however, is directly interpreted - without any inversion - as a response function 
(for more details see [3]). 



A. Simple example: deuteron photodisintegration 

In order to illustrate how the method works we first apply the LIT approach to a very 
simple physical problem, namely to the total deuteron photoabsorption cross section in 
unretarded dipole approximation. In this case the cross section is given by 

a^(a;) = in^auR'^iu) , (24) 

where a is the fine structure constant, uj is the energy of the photon absorbed by the 
deuteron, and R'^^{uj) denotes the response function defined as 

R» = ^KfimMoo - E., - E,) . (25) 

Here E^i and |0) are the deuteron bound state energy and wave function, while Enp and 
I/) denote relative kinetic energy and wave function of the outgoing np pair for a given 
two-nucleon Hamiltonian H: 

{H + E,)\0) = 0, {H - E^p)\f) = . (26) 

The transition operator B is defined by 

2 

where Zi and rf are the third components of position and isospin coordinates of the i-th 
nucleon. The LIT of is given by 

3 ^ 

L^((TK, cji) = = / R'iitu) Cico, aR, ai) du , (28) 

fc=i ^ 

where k = 1,2,3 correspond to different partial waves of the final state, namely ^Pq, ^Pi, 
and F2. 

First we consider deuteron photodisintegration with a realistic NN interaction. In the 



already mentioned review article of the LIT method IGj, such a case has been investigated 



using the AV14 NN potential 



251]. The L^ result is shown in Fig. 1, while the corresponding 
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FIG. 2: (T^ from inversion of of Fig. 1, up to 50 MeV (a) and in peak region (b), with various 
Mraax values: 10 (short dashed), 15 (dotted), 20 (long dashed), 25 (solid), 26 (dash-dotted). 

inversion results are illustrated in Fig. 2. For the inversion one observes a nice stability 
range of the results for all the shown Mmax values in the whole considered energy range, 
except for the peak region, where the inversion becomes stable only for higher Mmax- One 
notes that the Mmax values 25 and 26 lead essentially to identical results. 

In Fig. 3 the final inversion result (Mjnax=25) is compared with the corresponding cross 
section of a conventional calculation with explicit np continuum wave functions IGj. One 
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FIG. 3: Total deuteron photoabsorption cross section up to 50 MeV (a) and in peak region (b): 
LIT result (solid) and from calculation with explicit np continuum wave functions (crosses). 

finds an excellent agreement between the two calculations showing that one can reach high- 
precision results with the LIT method. 




E [MeVJ E [MeVJ 



FIG. 4: Phase shift ^Pi of fictitious np system at low (left) and higher (right) energies. 



Further LIT calculations for the deuteron total photoabsorption cross sections have al- 
ready been discussed in 26l . l27l| . For the aim of the present discussion, i.e. the way of 
working of the LIT method, 27| is particularly interesting. It is a case study for a fictitious 
np system with a low-lying and narrow resonance in the ^Pinucleon-nucleon partial wave 
(obtained by an additional attractive term, for details see 27|). The results of a conven- 
tional calculation with the fictitious np system for the ^Pi phase shifts and the "deuteron 
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FIG. 5: As Fig. 4 but for photodisintegration cross section. 

photoabsorption cross section" to the ^Pi final state are shown in Figs. 4 and 5, respectively. 
The ^Pi phase shifts exhibit two resonances, at Enp = 0.48 MeV and at about Enp = 10.5 
MeV. The low-energy resonance leads to the dominant structure of the photoabsorption 
cross section, a pronounced peak at a photon energy of 2.65 MeV with a width F of 270 
keV, while the second resonance only shows up as a rather tiny peak, which is more than 
four orders of magnitude smaller than the first peak. 
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FIG. 6: LIT L'i^{^Pi) of fictitious np system in first resonance region with aj = 1 MeV and various 
values of i?max (top) and ratio R=LJ^(i?max)/-^^[(-Rmax = 80 fm) (bottom). 

For the LIT calculation of the photoabsorption cross section of the fictitious np system 
the inversion basis set Xm (1151) is modified to account for the resonant structure, to this end 
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the functions Xm are relabeled: Xm — ^ Xm+i- In addition a new xi is defined, 

1/1 1 



{Enp - E,,,Y + (^)2 Vl + exp(-l) 1 + exp((E„p - a^)/a:i) 
where E^^s = ot^-, T = Q^s and 03 are additional nonlinear parameters. 




(0 [MeV] 
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FIG. 7: c7^(^Pi) of fictitious np system from inversion of L^|(^Pi, aj = 1 MeV, i?max) with various 
Afinax values and Rmax = 30 fm (top) and 50 fm (bottom) in first resonance region (left) and 
beyond (right); also shown results of a conventional calculation with explicit np continuum wave 
functions (crosses). 

The reason for such a LIT case study with a resonance lies in the results of a previous 
LIT calculation for the (e,e') longitudinal and transverse form factors of ^He {2^, where 
a resonance in the Coulomb monopole transition was obtained, but its width could not 
be determined. In the case study it is shown that for a proper resolution of a resonant 
structure it is very important to take into account the LIT function \^ up to rather large 
distances 0]. This has been checked (i) by solving ([6]) imposing at a two-nucleon distance 
r = -Rmax an asymptotic boundary condition which leads to a strong fall-off of \E'(^Pi) and 
(ii) by calculating the norm (\I'('^Pi)|\E'(^Pi)) only in the range from r = to r = -Rmax- In 
Fig. 6 we show the results for such a calculation choosing cr/ = 1 MeV. One notes that for a 
rather precise result, with errors below about 1%, one has to take -Rmax > 30 fm. A further 
increase of -Rmax to 50 fm leads to a reduction of the relative error by about a factor ten. In 
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Fig. 7 the inversions results with -Rmax = 30 and 50 fm are depicted in comparison to the 
result of the direct calculation. One observes that for both -Rmax values the resonance cross 
section is described with high accuracy. Differences between the two cases become evident 
in the region of the second maximum and at higher energy. In fact with -Rmax = 30 fm one 
finds only a reasonably good description, while a considerable improvement is obtained with 
-Rmax = 50 fm. 
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FIG. 8: As Fig. 6 but with cr/ = 0.5 MeV 

As next point we consider the reduction of aj from 1 MeV to 0.5 MeV. As shown in 
Fig. 8 one finds an enhancement of the relative error of by at least a factor of five in 
comparison to the corresponding -Rmax values of Fig. 6. The enhancement is easily under- 
stood by investigating the asymptotic solution of ([6]); in the here considered deuteron case 
it is described by the exponential fall-off exp{—r{Mai)^/h), where M is the nucleon mass. 
It is evident that a smaller aj leads to a longer range LIT function ^. As discussed in 
for cT/ = 0.1 MeV even -Rmax =300 fm is not completely sufficient, since only the resonance 
itself is resolved, but not the cross section at higher energies. 

The discussion above seems to infer that it is better to choose a rather large value for 
a J. On the other hand it should also be clear that the width of the resonance F and the 
value of aj are correlated. If aj is too large the resonance cannot be resolved. In the case 
study it has been found that cr/ = 2 MeV, about seven times larger than the width F, is 
still sufficient for a resolution of the resonance. However, in a general case the resonance 
width is not known beforehand and the question arises what should be the proper value for 
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(T/ in such a case. As pointed out in [27| one has to proceed as follows. One performs a LIT 
calculation with a given o"/ and compares the result to a LIT with a 5-peak in the response 
function or cross section. For example, in the here discussed deuteron case one sets 



The resulting (5-LIT is then given by the Lorentzian function 

L^(EP^^^ an, aj) = an, ai) , 

where Ri"^^^ is chosen such that the peak heights of L^ and L^J are equal. 
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FIG. 9: LIT L'^(^Pi) (full curves) and L^ (dashed curves) of fictitious np system in first resonance 
region with various cj/ values and -Rmax = 80 fm, except for aj = 0.1 MeV where i?max=300 fm. 

In order to obtain a reliable inversion the actual LIT L^J should have a larger width than 
the 5-LIT L^. If, on the contrary, they lead to essentially identical results, one has to reduce 
CT/ up to the point that the actual LIT is sufficiently different from the corresponding 5-LIT. 
In Fig. 9 we show such results for the deuteron case study. For cr/=5 and 10 MeV there are 
practically no differences between LIT and (5-LIT, while for o"/ < 2 MeV differences become 
visible. As a matter of fact aj=2 MeV is sufficient for a reliable inversion and thus one may 
conclude that in a general case one has to use a aj such that differences between LIT and 
5-LIT have at least the same size as in the aj=2 MeV case of Fig. 9. 
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In Fig. 10 we show the final inversion results from [271] with -Rmax =80 fm. One sees that 
the El resonance is precisely described with the LIT method for cr/ = 0.5 and 2 MeV, while 
the peak is somewhat underestimated with aj = 5 MeV. In the resonance region, with the 
two lower aj values, essentially the same results are obtained as for the cases of Fig. 7 with 
ai=l MeV and -Rmax =30 and 50 fm. From the comparison one further notes that the case 
(j/=2 MeV and -Rmax = 80 fm leads to even more precise results in the second resonance 
region, and beyond, than shown in Fig. 7 for -Rmax =50 fm. 




(0 [MeV] (0 [MeV] 



FIG. 10: As Fig. 7 but with -Rmax = 80 fm and ai = 5, 2, and 0.5 MeV. 



B. Reactions with A > 3 



Now we turn to realistic applications of the LIT method and consider first the ^He total 

h 

photoabsorption cross section. In fact this is one of the very first LIT applications [7[ , how- 
ever, initially only performed with semirealistic NN forces. In |3] it has been found that the 
calculated cross section shows a considerably more pronounced giant dipole resonance than 
the most recent experimental data of that time. Such a difference between experiment and 
an ab initio calculation has led to a renewed experimental interest in the ^He photoabsorp- 
tion cross section. In fact in the meantime three additional experiments have been carried 



out 



29 



30 



3l| . Finally, in js] the first calculation with realistic NN and 3N forces has also 
been published. In Fig. 11 we show the theoretical results of 5^ together with the experi- 
mental data. One notes that the 3N force leads to a considerable reduction of the peak cross 
section. One further sees that there is a very good agreement between theory and the data 
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of 3li and also quite a good agreement with the data of 29| , while the cross section of [30(] 
shows a completely different behavior. Also shown in the figure are data from experiments 
performed about 20 years ago. 
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IG. 11: Total ^He photoabsorption cross section: LIT calculation with AV18 NN potential only 
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It should be mentioned that today two other LIT calculations for the ^He total photoab- 



sorption cross section are available 
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where different realistic nuclear forces have been 



used. Essentially, they confirm the results of j^. 

In Fig. 12 LIT results for the ^Li and ^He total photoabsorption cross sections calculated 
with semirealistic NN forces are shown |^] . For ^Li one finds a single and rather broad cross 
section peak. On the contrary for ^He a very interesting cross section with a double peak 
becomes apparent in the calculation. This microscopic result can be interpreted as follows. 
In a cluster picture, with an a core and a di-neutron, the low-energy peak is due to the 
relative motion of di-neutron and a core. The second peak, however, cannot be obtained in 
a cluster model, but is explained by the classical El giant resonance picture with a collective 
response of all nucleons (relative motion of protons and neutrons). For ^Li, in a cluster model 
described by an a core plus a deuteron, a similar low-energy peak is missing, because the 
deuteron knock-out corresponds to an isoscalar transition, which cannot be induced by the 
isovector dipole operator (!27|) . On the other hand a transition to the antibound ^S'o(np) 
plus a core is possible. The nucleus ^Li exhibits a considerably larger width of the giant 
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dipole peak than ^He. In fact in the former a break-up into two three-body nuclei is possible 
(^H— ^He), while a similar break-up of ^He is not induced by the dipole operator, since the 
■^H— '^H pair has no dipole moment. The experimental ^He and ^Li photoabsorption cross 
sections are not yet well settled (see j^) and therefore not shown here. 
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FIG. 13: ^Li total photoabsorption cross section with semirealistic NN potential model AV4', 



experimental data from 



In Fig. 13 we depict the LIT calcu 
in comparison to experimental data 



ation for the ^Li total photoabsorption cross section 



42]. It is worthwhile to mention that the experimental 
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cross section has not been determined by summing up the various channel cross sections, 
but in a single experiment by the "diminution of photon flux" method. Like ^Li also ^Li has 
a giant dipole resonance peak with a rather large width. The comparison of experimental 
and theoretical results shows quite a good agreement, though only a semirealistic NN force 
has been used in the LIT calculation. Of course, it would be very interesting to have even 
more precise data and also a calculation with realistic nuclear forces. 




FIG. 14: RL{io,q) of (left) and ^He (right) at various q: AV18 NN potential (dotted), AV18 
NN + UIX 3N potential (solid); experimental data: triangles \^^, and circles [45]. 

Now we turn to the inclusive electron scattering response functions. For and ^He 
LIT calculations for the longitudinal response function Rrjco, q) have been carried out with 



43|, and above g=500 MeV/c 



realistic nuclear forces for momentum transfers below 
Relativistic corrections for the transition operator Ol have been taken into account and 
the frame dependence of the essentially nonrelativistic calculation has been studied. In 



431 ]. One notes that the 3N force reduces 



Fig. 14 Ri is shown for various lower q values 
the quasielastic peak height somewhat. The 3N force effect, however, does not lead to a 
consistent picture in comparison with experiment. In fact one flnds an improvement for ^He 
and a deterioration for '^H. 



19 




FIG. 15: RL{uJ,q) of ^He at various q, calculation with AV18 NN + UIX 3N potentials and two 



44| . triangles 



nucleon form factors fits: dipole (solid) and MMD from (dashed); experimental data: squares 



451 ] ■ and circles [4g]. 



In 



ll| it has been found that relativistic effects due to the kinetic energy can be largely 



reduced if Rl is first calculated in a specific reference frame, where the target nucleus moves 
with -Aq/2, and then transformed to the lab system. Different from a direct calculation in 
the lab system, one finds a correct result for the experimentally established quasielastic peak 
position ll|, as shown in Fig. 15. At g=500 and 600 MeV/c also for the peak height a good 



agreement with experimental data is obtained, whereas the theoretical Rl overestimates the 
experimental one at g=700 MeV/c (at even higher q experimental data are not available). 
As Fig. 15 shows also the choice of the nucleon form factor fit has a non-negligible impact 
on the result, but cannot explain the discrepancy with the data at g = 700 MeV/c. 

Recently we also calculated the transverse response function RT{uj,q) with realistic nu- 
clear forces {q < 500 MeV/c) 1^. Besides the usual one-body operators also tt and p 
exchange currents (EC), consistent with the NN potential model, were taken into account. 
In addition also the effect of the so-called Siegert operator has been studied. In Fig. 16 the 
theoretical results are shown together with experimental data. In the quasielastic region 
EC lead to some increase of the peak height (see left panels of Fig. 16), but the EC effect 
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FIG. 16: RT{uJ,q) in 



nRA NN potential 



quasielastic (left) and threshold region (right), theoretical results with Bon- 



481], and TM' 3N forces 



and following current operators: one-body (dotted) 



and one-body -|- vr-EC + p-EC -|- additional EC via Siegert operator (solid). Experimental data: 



left panels, triangles [44s], circles [45], and squares [46]; and right panels, circles 



is much larger close to threshold (see right panels of Fig. 16) and is important for a good 
agreement of theory and experiment. For g=500 MeV/c one finds different peak positions in 
theory and experiment, presumably due to the fact that Rt is calculated directly in the lab 
frame (see discussion above). The frame dependence oi Rt is presently under investigation. 

Presently we are extending our realistic calculation to the (e, e') response functions of 
^He and results for Rl will be published soon. 

Summarizing one may conclude that the LIT method is very powerful and allows a 
considerable extension of the range of microscopic ab initio calculations. 
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